version 6.4 



Charged and strange hadron elliptic flow 
in Cu+Cu collisions at v /s~^ 7 = 62.4 and 200 GeV 

B. I. Abelev, 8 M. M. Aggarwal, 31 Z. Ahammed, 48 A. V. Alakhverdyants, 18 I. Alekseev , 16 B. D. Anderson, 
D. Arkhipkin, 3 G. S. Averichev, 18 J. Balewski, 23 L. S. Barnby, 2 S. Baumgart, 53 D. R. Beavis, 3 R. Bellwied, 51 
J. Betancourt, 23 R. R. Bctts, 8 A. Bhasin, 17 A. K. Bhati, 31 H. Bichsel, 50 J. Bielcik, 10 J. Bielcikova, 11 B. Biritz 

L. C. Bland, 3 B. E. Bonner, 37 J. Bouchct, 19 E. Braidot, 28 A. V. Brandin, 26 A. Bridgeman, 1 E. Bruna, 53 
S. Bueltmann, 30 I. Bunzarov, 18 T. P. Burton, 3 X. Z. Cai, 41 H. Caines, 53 M. Calderon de la Barca Sanchez, 5 
O. Catu, 53 D. Cebra, 5 R. Cendejas, 6 M. C. Cervantes, 43 Z. Chajecki, 29 P. Chaloupka, 11 S. Chattopadhyay, 41 
H. F. Chen, 39 J. H. Chen, 41 J. Y. Chen, 52 J. Cheng, 45 M. Cherney, 9 A. Chikanian, 53 K. E. Choi, 35 W. Christi 
P. Chung, 11 R. F. Clarke, 43 M. J. M. Codrington, 43 R. Corliss, 23 J. G. Cramer, 50 H. J. Crawford, 4 
D. Das, 5 S. Dash, 13 A. Davila Leyva, 44 L. C. Dc Silva, 51 R. R. Debbe, 3 T. G. Dedovich, 18 M. DePhillips, 3 
A. A. Derevschikov, 33 R. Derradi de Souza, 7 L. Didcnko, 3 P. Djawotho, 43 S. M. Dogra, 17 X. Dong, 22 
J. L. Drachcnberg, 43 J. E. Draper, 5 J. C. Dunlop, 3 M. R. Dutta Mazumdar, 48 L. G. Efimov, 18 E. Elhalhuli, 2 
M. Elnimr, 51 J. Engelage, 4 G. Eppley, 37 B. Erazmus, 42 M. Estiennc, 42 L. Eun, 32 O. Evdokimov, 8 P. Fachini, 3 
R. Fatcmi, 20 J. Fedorisin, 18 R. G. Fersch, 20 P. Filip, 18 E. Finch, 53 V. Fine, 3 Y. Fisyak, 3 C. A. Gagliardi, 43 
D. R. Gangadharan, 6 M. S. Ganti, 48 E. J. Garcia-Solis, 8 A. Geromitsos, 42 F. Gcurts, 37 V. Ghazikhanian, 6 
P. Ghosh, 48 Y. N. Gorbunov, 9 A. Gordon, 3 O. Grebenyuk, 22 D. Grosnick, 47 B. Grube, 35 S. M. Guertin, 6 

A. Gupta, 17 N. Gupta, 17 W. Guryn, 3 B. Haag, 5 A. Hamed 43 L-X. Han, 41 J. W. Harris, 53 J. P. Hays-Wehle, 23 
M. Heinz, 53 S. Heppelmann, 32 A. Hirsch, 34 E. Hjort, 22 A. M. Hoffman, 23 G. W. Hoffmann, 44 D. J. Hofman, 8 
R. S. Hollis, 8 H. Z. Huang, 6 T. J. Humanic, 29 L. Huo, 43 G. Igo, 6 A. Iordanova, 8 P. Jacobs, 22 W. W. Jacobs, 15 

P. Jakl, 11 C. Jena, 13 F. Jin, 41 C. L. Jones, 23 P. G. Jones, 2 J. Joseph, 19 E. G. Judd, 4 S. Kabana, 42 K. Kajimoto, 44 
K. Kang, 45 J. Kapitan, 11 K. Kauder, 8 D. Keane, 19 A. Kechechyan, 18 D. Kettler, 50 D. P. Kikola, 22 J. Kiryluk, 22 
A. Kisiel, 49 S. R. Klein, 22 A. G. Knospe, 53 A. Kocoloski, 23 D. D. Koctke, 47 T. Kollegger, 12 J. Konzer, 34 
M. Kopytine, 19 I. Koralt, 30 L. Koroleva, 16 W. Korsch, 20 L. Kotchenda, 26 V. Kouchpil, 11 P. Kravtsov, 26 
K. Krueger, 1 M. Krus, 10 L. Kumar, 31 P. Kurnadi, 6 M. A. C. Lamont, 3 J. M. Landgraf, 3 S. LaPointe, 51 J. Lauret, 3 

A. Lebedev, 3 R. Lednicky, 18 C-H. Lee, 35 J. H. Lee, 3 W. Leight, 23 M. J. LeVine, 3 C. Li, 39 L. Li, 44 N. Li, 52 
W. Li, 41 X. Li, 40 X. Li, 34 Y. Li, 45 Z. Li, 52 G. Lin, 53 S. J. Lindenbaum, 27 'E M. A. Lisa, 29 F. Liu, 52 H. Liu, 5 

J. Liu, 37 T. Ljubicic, 3 W. J. Llopc, 37 R. S. Longacre, 3 W. A. Love, 3 Y. Lu, 39 G. L. Ma, 41 Y. G. Ma, 41 

D. P. Mahapatra, 13 R. Majka, 53 O. I. Mall, 5 L. K. Mangotra, 17 R. Manweiler, 47 S. Margetis, 19 C. Markert, 44 
H. Masui, 22 H. S. Matis, 22 Yu. A. Matulenko, 33 D. McDonald, 37 T. S. McShane, 9 A. Meschanin, 33 R. Milner, 23 

N. G. Minaev, 33 S. Mioduszewski, 43 A. Mischke, 28 M. K. Mitrovski, 12 B. Mohanty, 48 M. M. Mondal, 48 

B. Morozov, 16 D. A. Morozov, 33 M. G. Munhoz, 38 B. K. Nandi, 14 C. Nattrass, 53 T. K. Nayak, 48 J. M. Nelson, 2 
P. K. Netrakanti, 34 M. J. Ng, 4 L. V. Nogach, 33 S. B. Nurushcv, 33 G. Odynicc, 22 A. Ogawa, 3 H. Okada, 3 

V. Okorokov, 26 D. Olson, 22 M. Pachr, 10 B. S. Page, 15 S. K. Pal, 48 Y. Pandit, 19 Y. Panebratsev, 18 T. Pawlak 49 
T. Peitzmann, 28 V. Percvoztchikov, 3 C. Perkins, 4 W. Peryt, 49 S. C. Phatak, 13 P. Pile, 3 M. Planinic, 54 
M. A. Ploskon, 22 J. Pluta, 49 D. Plyku, 30 N. Poljak, 54 A. M. Poskanzcr, 22 B. V. K. S. Potukuchi, 17 C. B. Powell, 22 
D. Prindlc, 50 C. Pruneau, 51 N. K. Pruthi, 31 P. R. Pujahari, 14 J. Putschke, 53 R. Raniwala, 36 S. Raniwala, 36 
R. L. Ray, 44 R. Rcdwinc, 23 R. Reed, 5 J. M. Rehberg, 12 H. G. Ritter, 22 J. B. Roberts, 37 O. V. Rogachcvskiy, 18 

J. L. Romero, 5 A. Rose, 22 C. Roy, 42 L. Ruan, 3 R. Sahoo, 42 S. Sakai, 6 I. Sakrcjda, 22 T. Sakuma, 23 S. Salur, 5 
J. Sandweiss, 53 E. Sangaline, 5 J. Schambach, 44 R. P. Scharenberg, 34 N. Schmitz, 24 T. R. Schuster/ 2 J. Seele, 23 
J. Seger, 9 I. Selyuzhenkov, 15 P. Seyboth, 24 E. Shahaliev, 18 M. Shao, 39 M. Sharma, 51 S. S. Shi, 52 >0X. H. Shi, 41 

E. P. Sichtermann, 22 F. Simon, 24 R. N. Singaraju, 48 M. J. Skoby, 34 N. Smirnov, 53 P. Sorensen, 3 J. Sowinski, 15 
H. M. Spinka, 1 B. Srivastava, 34 T. D. S. Stanislaus, 47 D. Staszak, 6 J. R. Stevens, 15 R. Stock, 12 M. Strikhanov, 26 

B. Stringfellow, 34 A. A. P. Suaide, 38 M. C. Suarez, 8 N. L. Subba, 19 M. Sumbera, 11 X. M. Sun, 22 Y. Sun, 39 
Z. Sun, 21 B. Surrow, 23 D. N. Svirida, 16 T. J. M. Symons, 22 A. Szanto de Toledo, 38 J. Takahashi, 7 A. H. Tang, 3 

Z. Tang, 39 L. H. Tarini, 51 T. Tarnowsky, 25 D. Thcin, 44 J. H. Thomas, 22 J. Tian, 41 A. R. Timmins, 51 
S. Timoshenko, 26 D. Tlusty, 11 M. Tokarev, 18 V. N. Tram, 22 S. Trcntalange, 6 R. E. Tribble, 43 O. D. Tsai, 6 
J. Ulery, 34 T. Ullrich, 3 D. G. Underwood, 1 G. Van Buren, 3 M. van Leeuwen, 28 G. van Nieuwenhuizen, 23 
J. A. Vanfosscn, Jr., 19 R. Varma, 14 G. M. S. Vasconcelos, 7 A. N. Vasiliev, 33 F. Videback, 3 Y. P. Viyogi 48 
S. Vokal, 18 S. A. Voloshin, 51 M. Wada, 44 M. Walker, 23 F. Wang, 34 G. Wang, 6 H. Wang, 25 J. S. Wang, 21 
Q. Wang, 34 X. L. Wang, 39 Y. Wang, 45 G. Webb, 20 J. C. Webb, 3 G. D. Westfall, 25 C. Whittcn Jr., 6 



5 



;ie 3 



2 



H. Wieman, 22 E. Wingficld, 44 S. W. Wissink, 15 R. Witt, 46 Y. Wu, 52 W. Xie, 34 N. Xu, 22 Q. H. Xu 40 
W. Xu, 6 Y. Xu, 39 Z. Xu, 3 L. Xue, 41 Y. Yang, 21 P. Yepcs, 37 K. Yip, 3 I-K. Yoo, 35 Q. Yuc 45 M. Zawisza, 49 
H. Zbroszczyk 49 W. Zhan, 21 S. Zhang, 41 W. M. Zhang, 19 X. P. Zhang, 22 Y. Zhang, 22 Z. P. Zhang, 39 J. Zhao, 41 
C. Zhong 41 J. Zhou, 37 W. Zhou, 40 X. Zhu 45 Y. H. Zhu 41 R. Zoulkarneev, 18 and Y. Zoulkarneeva 18 

(STAR Collaboration) 

1 Argonne National Laboratory, Argonne, Illinois 60439, USA 
2 University of Birmingham, Birmingham, United Kingdom 
3 Brookhaven National Laboratory, Upton, New York 11973, USA 
4 University of California, Berkeley, California 94720, USA 
5 University of California, Davis, California 95616, USA 
6 University of California, Los Angeles, California 90095, USA 

7 Universidade Estadual de Campinas, Sao Paulo, Brazil 
8 University of Illinois at Chicago, Chicago, Illinois 60607, USA 
9 Creighton University, Omaha, Nebraska 68178, USA 
10 Czech Technical University in Prague, FNSPE, Prague, 115 19, Czech Republic 
11 Nuclear Physics Institute AS CR, 250 68 Rez/Prague, Czech Republic 
University of Frankfurt, Frankfurt, Germany 
1 3 Institute of Physics, Bhubaneswar 751005, India 
14 Indian Institute of Technology, Mumbai, India 
15 Indiana University, Bloommgton, Indiana 47408, USA 
1 ' Alikhanov Institute for Theoretical and Experimental Physics, Moscow, Russia 
17 University of Jammu, Jammu 180001, India 
18 Joint Institute for Nuclear Research, Dubna, 141 980, Russia 
"'Kent State University, Kent, Ohio 44242, USA 
20 University of Kentucky, Lexington, Kentucky, 40506-0055, USA 
21 Institute of Modern Physics, Lanzhou, China 
22 Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 
23 Massachusetts Institute of Technology, Cambridge, MA 02139-4307, USA 
24 Max- Planck- Institut fur Physik, Munich, Germany 
25 Michigan State University, East Lansing, Michigan 48824, USA 
26 Moscow Engineering Physics Institute, Moscow Russia 
27 City College of New York, New York City, New York 10031, USA 
28 NIKHEF and Utrecht University, Amsterdam, The Netherlands 
20 Ohio State University, Columbus, Ohio 43210, USA 
30 Old Dominion University, Norfolk, VA, 23529, USA 
31 Panjab University, Chandigarh 160014, India 
32 Pennsylvania State University, University Park, Pennsylvania 16802, USA 
Institute of High Energy Physics, Protvino, Russia 
34 Purdue University, West Lafayette, Indiana 47907, USA 
S5 Pusan National University, Pusan, Republic of Korea 
36 University of Rajasthan, Jaipur 302004, India 
37 Rice University, Houston, Texas 77251, USA 
38 Universidade de Sao Paulo, Sao Paulo, Brazil 
39 University of Science & Technology of China, Hefei 230026, China 
4 °Shandong University, Jinan, Shandong 250100, China 
4 Shanghai Institute of Applied Physics, Shanghai 201800, China 
42 SUBATECH, Nantes, France 
43 Texas A&M University, College Station, Texas 77843, USA 
44 University of Texas, Austin, Texas 78712, USA 
45 Tsmghua University, Beijing 100084, China 
46 United States Naval Academy, Annapolis, MD 21402, USA 

47 Valparaiso University, Valparaiso, Indiana 46383, USA 
48 Variable Energy Cyclotron Centre, Kolkata 700064, India 
49 Warsaw University of Technology, Warsaw, Poland 
50 University of Washington, Seattle, Washington 98195, USA 
51 Wayne State University, Detroit, Michigan 48201, USA 
52 Institute of Particle Physics, CCNU (HZNU), Wuhan 430079, China 
53 Yale University, New Haven, Connecticut 06520, USA 
54 University of Zagreb, Zagreb, HR-10002, Croatia 
(Dated: March 26, 2010) 



3 



We present the results of an elliptic flow, V2, analysis of Cu+Cu collisions recorded with the STAR 
detector at RHIC at ^/snn — 62.4 and 200 GeV. Elliptic flow as a function of transverse momen- 
tum, V2(pt), is reported for different collision centralities for charged hadrons h , and strangeness 
containing hadrons Kg, A, S, (j> in the midrapidity region \rj\ < 1.0. Significant reduction in sys- 
tematic uncertainty of the measurement due to non-flow effects has been achieved by correlating 
particles at midrapidity, \rj\ < 1.0, with those at forward rapidity, 2.5 < \r]\ < 4.0. We also present 
azimuthal correlations in p+p collisions at i/s = 200 GeV to help estimating non-flow effects. To 
study the system-size dependence of elliptic flow, we present a detailed comparison with previously 
published results from Au+Au collisions at ^/snn = 200 GeV. We observe that vzipr) of strange 
hadrons has similar scaling properties as were first observed in Au+Au collisions, i.e.: (i) at low 
transverse momenta, pr < 2 GeV /c, V2 scales with transverse kinetic energy, mr — m, and (ii) at 
intermediate pr, 2 < pr < 4 GeV/c, it scales with the number of constituent quarks, n q . We have 
found that ideal hydrodynamic calculations fail to reproduce the centrality dependence of V2(pr) 
for K s and A. Eccentricity scaled V2 values, V2/e, are larger in more central collisions, suggesting 
stronger collective flow develops in more central collisions. The comparison with Au+Au collisions 
which go further in density shows V2/E depend on the system size, number of participants N pa , It . 
This indicates that the ideal hydrodynamic limit is not reached in Cu+Cu collisions, presumably 
because the assumption of thermalization is not attained. 

PACS numbers: 25.75.Ld, 25.75.Dw 



I. INTRODUCTION 



At the early stages of rclativistic heavy ion collisions, 
a hot and dense, strongly interacting medium is created. 
The subsequent system evolution is determined by the 
nature of the medium. Experimentally, the dynamics 
of the system evolution has been studied by measuring 
the azimuthal anisotropy of the particle production rel- 
ative to the reaction plane [J-Q- The centrality of the 
collision, defined by the transverse distance between the 
centers of the colliding nuclei called the impact parame- 
ter, results in an "almond-shaped" overlap region that 
is spatially azimuthal anisotropic. It is generally as- 
sumed that the initial spatial anisotropy in the system 
is converted into momentum-space anisotropy through 
re-scatterings [J, [f|. The elliptic flow, u 2 , is the second 
harmonic coefficient of a Fourier expansion of the final 
momentum-space azimuthal anisotropy. Due to the self- 
quenching effect, it provides information about the dy- 
namics at the early stage of the collisions . Elliptic 
flow can provide information about the pressure gradi- 
ents, the effective degrees of freedom, the degree of ther- 
malization, and equation of state of the matter created 
at the early stage. Thus, the centrality and system-size 
dependence of elliptic flow at different beam energies can 
be used to study theproperties of the matter created in 
heavy ion collisions [f| . 

Recently, two important insights have been obtained 
from the experimental results on V2 as a function of 
transverse momentum, pr, in Au+Au collisions at RHIC. 
First, in the low px region, px < 2 GeV/c, the hadron 
mass hierarchy predicted by ideal hydrodynamic calcu- 



lations is observed for identified hadrons 7r, K, Kg, p, 
A and 5 0-[l|. Even the 4> an d f2, which are be- 
lieved to have a reduced cross section for hadronic in- 
teractions I14M19I1 . are consistent with the mass order- 
ing [ID, [20L |2l|. Second, in the intermediate pr region, 
2 < pt < 4 GeV/c, i>2 (pt) follows a scaling depend- 
ing on the number of constituent quarks within a given 
hadron, which can be explained via coalescence mod- 
els [HI, I20I [22j . Quark number scaling suggests that the 
system is in a partonic state and that the constituent 
quark degrees of freedom were relevant during the time 
V2 was developed. 

STAR's first published paper showed that elliptic flow 
at RHIC is unexpectedly large [Hf , comparable to p redic- 
tions of ideal hydrodynamic calculations [3. |24| - 126{ . This 
observation is among the evidence favoring the picture of 
a nearly-perfect liquid [13] • With the assumption of ther- 
malization, ideal hydrodynamic calculations predict that 
the V2 divided by spatial eccentricity, e, does not depend 
on the collision centrality [28|. The spatial eccentricity 
is denned by [H]: 



Or 



W 2 



(i) 
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where x and y are the spatial coordinates in the plane 
perpendicular to the collision axis. The angle brackets ( ) 
denote an average weighted by the initial density. How- 
ever, recent RHIC v^js data for charged hadrons hr- and 
strangeness containing hadrons Kg, <fi, A, S show a trend 
to increase as a function of the particle density scaled by 
the system-size [3(J HH , lacking the saturation indicated 
by ideal hydrodynamic calculations [3l| . This monotonic 
increase is a feature of a class of model descriptions that 
conform to the low-density limit [H| ■ Whether the ther- 
malization and ideal hydrodynamic limit are reached or 
not at RHIC is not conclusive. A transport model sug- 
gested in Ref. [H| is constructed to link the low-density 
limit to the ideal hydrodynamic limit. In the micro- 
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scopic transport picture, the ideal hydrodynamic limit 
is reached when the mean free path is very small or the 
cross section is very large. With this transport model 
approach, the degree of thermalization and the ideal hy- 
drodynamic limit can be addressed J34[ ■ 

The previous results mainly focus on the centrality de- 
pendence of charged hadron and identified hadron vi in 
Au+Au collisions. Since the conditions in Au+Au col- 
lisions might not hold in smaller systems and at lower 
beam energies, the system-size and beam-energy depen- 
dence of identified hadron elliptic flow will shed light on 
the properties of partonic collectivity and quark degrees 
of freedom. Further, the study of i>2 in collisions of nu- 
clei smaller than Au+Au will allow us to test the early 
thermalization hypothesis in Au+Au collisions. To date, 
there are only a few studies of identified hadron «2 in 
Cu+Cu collisions. In this article, we present the results 
on the azimuthal anisotropy parameter V2(pt) of h^ 1 , Kg, 
A , S and (j) from ^fs^ = 62.4 and 200 GeV Cu+Cu 
collisions. As a function of collision centrality, the scal- 
ing properties of V2 with the transverse kinetic energy 
rax — m and the number of constituent quarks n q are re- 
ported. In the quantity = \Jp\ + m 2 , m denotes the 
rest mass of a given hadron. We also discuss system-size 
dependence in this article. 

The rest of the paper is organized in the following 
way: Section |n] summarizes the analysis details includ- 
ing data and centrality selections, particle identification 
and flow methods used for charged hadrons and identi- 
fied hadrons. In the following, we use h , A and 3 to 
denote charged hadron, A + A and S~ +S , respectively. 
In Section IlIIl we present measurements of vi for h in 
Cu+Cu collisions from different analysis methods. Dif- 
ferences in V2 are used to estimate the systematic error. 
Section [TV] presents the results and physics discussion of 
the scaling properties and system-size dependence along 
with ideal hydrodynamic calculations. Last, a summary 
is presented in Section [V] 



II. EXPERIMENTS AND ANALYSIS 
A. Experiments and data sets 

For this article, our data were collected from ^/snn 
= 62.4 and 200 GeV Cu+Cu collisions with the STAR 
detector during the fifth RHIC run in 2005. In addition 
data from ^/snn = 200 GeV p+p collisions in 2005 were 
used in the analysis of non-flow contributions. STAR's 
main Time Projection Chamber (TPC) (HJ and two For- 
ward Time Projection Chambers (FTPCs) [36[ were used 
for particle tracking in the central region (|?7| < 1.0) and 
forward regions (2.5 < |?7| < 4.0), respectively. Both the 
TPC and FTPCs provide azimuthal acceptance over 27r. 
Only those events which have the primary vertex position 
along the longitudinal beam direction (V z ) within 30 cm 
of the nominal collision point are selected for the analy- 
sis. This is done in order to have a more uniform detector 



performance within |7/| < 1.0. The centrality definition, 
which is based on the raw charged particle TPC multi- 
plicity with | ?y | < 0.5, is the same as used previously [37| . 
After quality cuts, the number of the 60% most central 
events is about 24 million for 200 GeV Cu+Cu collisions 
and 10 million for 62.4 GeV Cu+Cu collisions. The re- 
sults from more peripheral collisions are not presented 
due to trigger inefficiencies at low multiplicity. 

The centrality was defined using the number of charged 
tracks with quality cuts similar to those in Ref. [U . The 
60% most central events for i> 2 analysis of were di- 
vided into six centrality bins, each spanning an interval 
of 10% of the geometric cross section. For V2 analysis 
of K s and A, centrality bins of - 20% and 20 - 60% 
were used. In order to reduce the multiplicity fluctua- 
tions in wide centrality bins, we calculated V2 in the 10% 
wide bins, then combined them using particle yield as the 
weight. 

To select good tracks from primary collisions, charged 
particle tracks coming from the collision which trans- 
versed the TPC or FTPCs were selected by requiring the 
distance of closest approach to the primary vertex be less 
than 3 cm. Tracks used for K s , A and H reconstruction 
were not subject to this cut. We required that the TPC 
and FTPCs had a number of hits used for reconstruction 
of the tracks of the particles > 15 and > 5, respectively. 
For the TPC and FTPCs the ratio of the number of fit 
hits to maximum possible hits was > 0.52. An additional 
transverse momentum cut (0.15 < pr < 2 GeV/c) was 
applied to the charged tracks for the event plane deter- 
mination. 



B. Particle identification 

We utilized the topology of decay as measured with 
the TPC inside the magnetic field to identify Kg, <j>, 



A (A) and 5 (^ + )- We used the decay channels: 
K s -> 7r++7r~, <j> -> K + +K~, A -> p+Tf- (A -> P+7T+), 

S~ — > A + 7r~ (S + — > A + 7r + ). Similar to the previous 
analysis in Ref. [3l| , topological and kinematic cuts were 
applied to reduce the combinatorial backgrounds. Fig- 
ure [T] shows the invariant mass distributions for (a) Kg, 
(b) (p, (c) A and (d) S for selected pr bins in y/s^N 
= 200 GeV Cu+Cu 60% most central collisions. The 
cuts used for Cu+Cu collisions are similar to those for 
Au+Au collisions in Ref. (3l|. The combinatorial back- 
grounds were estimated from the fourth order polynomial 
The invariant mass distribution 
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fits for Kg 1 and A 
for <j) is after subtraction of combinatorial background 
estimated using event mixing [38|, the remaining com- 
binatorial backgrounds were estimated by a first order 
polynomial fit [20(. For 5, the background was estimated 
by rotating the transverse momentum of the daughter A 
by 180°. This operation breaks the correlation between 
the A and the other daughter particle. The resulting 
invariant mass distributions provide a good approxima- 
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Invariant Mass (GeV/c 2 ) 



FIG. 1: (Color online) Invariant mass distributions for (a) K° s (1.2 < p T < 1.4 GeV/c), (b) cf> (1.0 < p T < 2.0 GeV/c), (c) 
A (1.4 < pt < 1.6 GeV/c) and (d) 3 (1.25 < p T < 1.75 GeV/c) in y4^7 = 200 GeV Cu+Cu 60% most central collisions. 
The solid curves represent the fits to the invariant mass distributions: Gaussians plus fourth-order polynomials for K s , A and 
3, Breit-Wigner plus a linear function for cj>- The dotted curves are the estimated backgrounds: the fourth order polynomials 
for Kg and A, a linear function for cj>, and a rotation method described in the text for 3. For clarity, the invariant mass 
distributions for K° s , A, (j> and 3 are scaled by 1/50 000, 1/130 000, 1/5 000 and 1/8 000, respectively. The error bars are 
shown only for the statistical uncertainties. 



tion of the true background distribution. The detailed 
description of the method can be found in Refs. (l2l. [l3|. 



C. Flow methods 

Anisotropic flow, which is an anisotropy in the parti- 
cle production relative to the reaction plane, results in 
correlations among particles and can be studied by the 
analysis of these correlations. At the same time these 
correlations are affected by other effects which are not 
related to the orientation of the reaction plane. Such are 
commonly referred to as non-flow, and are due, e.g., to 
resonance decays and jet production. Different methods 
used to measure anisotropic flow are affected by non-flow 
effects in different ways, and are used in this analysis to 
evaluate the systematic uncertainty of the measurements. 



1. Event Plane method with TPC event plane 

The Event Plane method Q uses the anisotropic flow 
itself to determine the event plane (the estimated reac- 
tion plane), which can be done for each harmonic. The 
second harmonic flow vector, Q 2 , of the event is con- 
structed using the TPC tracks i in the event with their 
azimuthal angle, fa, according to Eq. ^ and Eq. ([3]). In 
order to maximize the resolution of the flow effect, the 
weights Wi are set equal to pt up to 2 GeV /c. 

Q 2 cos(2^ 2 ) = Qix = ^2wiCos(2fa) (2) 



Q 2 sin(2* 2 ) = Q2 V = ^2wism(2fa) (3) 

Elliptic flow is first calculated with respect to the event 
plane angle $2 a s shown in Eq. which is called the 
observed w 2 . The angle brackets indicate an average over 
all particles in all events. However, tracks used for the 
V2 calculation arc excluded from the calculation of the 
flow vector to remove auto-correlation effects. Then the 
observed i> 2 is corrected by the event plane resolution 
(the denominator in Eq. ([5)1) to obtain w 2 relative to the 
event plane. 

v? s = (cos[2(0-* 2 )]) (4) 



V ' 2 (C0S[2(* 2 - ttr)]) 1 ' 

The results are denoted as u 2 {TPC} in the following. 

Since the reaction plane is unknown, the denominator 
in Eq. ([5]) is still not calculable. As shown in Eq. ©, we 
estimate the event plane resolution by the correlations 
between the azimuthal angles of two subset groups of 
tracks, called subevents A and B. In this analysis, we use 
two random subevents with equal numbers of particles. 
In Eq. J5| C is a constant calculated from the known 
multiplicity dependence of the resolution 0). 

(cos[2(* a -* r )]) = C^/(cos[2(^ - *f )]) (6) 

In the case of low resolution (< 0.2), such as for the 
FTPC event plane, C approaches s/2. 
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The reaction plane azimuthal distribution should be 
isotropic in the laboratory frame. Thus, the event plane 
azimuthal distribution has to be flat if the detectors have 
ideal acceptance. Since the detectors usually have non- 
uniform acceptance, a procedure for flattening the event 
plane distribution is necessary. For the event plane re- 
constructed from TPC tracks, the (ft weight method is an 
effective way to flatten the distribution. (ft weights are 
generated by inverting the (ft distributions of detected 
tracks for a large event sample. The detector accep- 
tance bias is removed by applying the (ft weight at the 
(ft of each track to that track. The (ft weights are folded 
into the weight Wi in Eq. ([2]) and Eq. ([3]). Independent 
corrections were applied to each centrality selection in 
10% increments and in 2 bins in the primary vertex po- 
sition along the longitudinal beam direction (V z ). The 
corrections were done on a run by run basis (around 50 k 
events). 



2. Event Plane method with FTPC event plane 

The rj gap between two FTPCs sitting at two sides of 
the collision in the forward regions can be used to reduce 
non-flow effects due to short-range correlations. The ba- 
sic procedures are similar to those for the Event Plane 
method with the TPC event plane. There are three steps: 
estimate the event plane with FTPC tracks, calculate V2 
with respect to the event plane, and obtain the real vi by 
correction to the real reaction plane. Equations ©-© 
can be applied, except that: i) the sums in Eq. and 
Eq. © go over FTPC tracks instead of TPC tracks, and 
ii) two subset groups of tracks are classified according 
to the sign of rj. The tracks with —4 < 77 < —2.5 and 
2.5 < rj < 4 are called East subevent and West subevent, 
respectively. Hence, the resolution in Eq. ([6]) is calcu- 
lated by the correlation between the azimuthal angles 
*f ast and \l/ Wost Tne average m E q gj runs over the 

TPC tracks as before. The result of this procedure is 
denoted as w 2 {FTPC}. 

Due to the serious loss of acceptance for FTPCs due to 
partially non-functioning readout electronics, the num- 
ber of tracks detected by the best sector is about 6 
times greater than for the worst one. The result is that 
the <j> weight method is not enough to generate a flat 
event plane distribution. Thus, further small corrections 
are applied after (ft weight corrections using the shift 
method j3{|. Equation ([7]) shows the formula for the 
shift correction. The averages in Eq. ([7]) are taken from 
a large sample of events. In this analysis, the correction 
is done up to the twentieth harmonic. This was done in 
order to make the \ 2 divided by the number of degrees 
of freedom of a flat fit to the event plane azimuthal an- 
gle distribution to be less than 1. The distributions of 
\E»f ast and *^ ost are separately flattened and then the 
full-event event plane distribution is flattened. Accord- 
ingly, the observed V2 and resolution are calculated using 
the shifted (sub)cvent plane azimuthal angles. 



= * + £ n ±[-(sin(2ratf))cos(2ntf) 
+ <cos(2n*)) sin(2n*)] 



(7) 



3. Scalar Product method 

The Scalar Product method [3(1 S3] is similar to the 
Event Plane method, and gives i>2 as: 



v 2 (pt) 



(Q 2 «2,i(Pr)) 



(8) 



where U2.i = cos(2cfti) + ism(2(fti) is a unit vector of the 
ith particle, Q2 = Ylk * s ^ ne ^ ow vec tor with the 
sum running over all other particles k in the event. The 
superscript * denotes the complex conjugate of a complex 
number. A and B denote the two subevents. In the case 
that Q2 is normalized to a unit vector, Eq. (|8]) reduces to 
the Event Plane method. In the Scalar Product method, 
one can use a different (re-centering) technique [Il| to 
correct for detector effects, which presents an alterna- 
tive to the weighting and shifting procedures described 
in Sec. Ill C 11 2 above. The Scalar Product method is 
applied to the V2 measurement of charged hadrons. 



D. method 

For V2 of the identified particles K% , (ft, A and S, the 
i>2 versus TOi nv method is used [111, |42| . Since V2 is addi- 
tive, one can write the total i> 2 lg+Bs as a sum of Signal 
and Background contributions weighted by their relative 
yields: 



v 



Sig+Bg 



(TOinv) 



Sig 



Sig 



' (TOinv) 



Sig + Bg 
Bg 



Sig + Bg 



(TOinv ) + 
(TOinv) 



(9) 



This method involves the calculation of v 2 lg+B& as a func- 
tion of TOinv and then fitting the distribution using Eq. (|9]) 
with measured relative yields and parameterizations of 
i>2 Ig and «2 S ' The sig+Bg ( m fav) distribution is the Bg 
divided by (Sig + Bg). The Si g^ s Bo - (TO inv ) distribution is 

simply calculated by 1 — gig+'Bg ( TO inv)- The term is 
parameterized as a linear function in order to take care 
of the non-constant value as a function of mi nv . The 
fit result w 2 Ig is the final observed v-z. Why this method 
works well for measuring signal i>2 is explained as follows: 
a set of data points is used in the fit over a wide mi nv 
region for Sig and Bg. Data points far from the mass 
peak constrain i^^mhiv), since pure Bg is expected in 
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FIG. 2: (Color online) Charged hadron azimuthal correlations 
as a function of pr in ^snn = 200 GeV 60% most central 
Cu+Cu collisions (closed squares) compared to those from 
t/snn = 200 GeV p+p collisions (open squares) . Flow vector 
calculated from (a) TPC tracks, (b) FTPC tracks. The error 
bars are shown only for the statistical uncertainties. 



this region. Under the peak, the v. 



Sig+Bg 



(m- mv ) is dom- 



inated by the Sig distribution. Finally, the v 2 signal is 
extracted by the fitting method shown in Eq. (|9]). 

Note that the subtraction procedure used to extract 
the Vi signal for a given identified particle is indepen- 
dent of the flow correlations. The i> 2 distributions of the 
overall signal and background are evaluated by one of the 
flow analysis methods discussed in Sec. Ill C 11 3. In this 
paper, the Event Plane method with the FTPC event 
plane is applied for Kg, (j>, A and S. 



E. Non-flow contribution for various methods 

The method of determining v 2 using cumulants of var- 
ious orders has been shown to eliminate non-flow cor- 
relations. However, the method is useful only for large 
values of flow and multiplicity. For the relatively low val- 
ues of flow and multiplicity seen in Cu+Cu collision, the 
non-flow correlations have been estimated, as described 
below. 

The Event Plane method with the TPC event plane is 
sensitive to non-flow effects. Particles of interest tend to 
correlate with particles used in the flow vector calculation 
due to short-range non-flow correlations. Also, particles 
of two random sub-events tend to have those correlations. 
Thus, non-flow exists in both the observed v 2 (Eq. ^) 



FIG. 3: (Color online) (a) Charged hadron V2(pt) in ^/snn 
— 200 GeV — 60% Cu+Cu collisions. Open circles, closed 
circles, open squares and closed squares represent the results 
of V2 as function of px measured by the TPC flow vector 
(« 2 {TPC}), the FTPC flow vector (u 2 {FTPC}), the TPC and 
FTPC flow vector with subtracting the azimuthal correlations 
in p+p collisions (v 2 {AA - pp, TPC}, V2{AA - pp, FTPC}). 
(b) The ratio of the results for the various methods described 
in (a). The error bars are shown only for the statistical un- 
certainties. 



and the resolution (Eq. ([5])). To reduce non-flow effects 
due to short-range correlations, we take advantage of the 
large 77 gap between the two FTPCs sitting at the two 
sides of the collision in the forward regions. Non-flow 
is reduced by the r\ gap between the TPC and FTPCs, 
but this may not be large enough to remove all non-flow 
correlations. Thus, we investigate these effects by com- 
paring the azimuthal correlations measured in Cu+Cu 
to those in p+p collisions, where all correlations are as- 
sumed to be of non-flow origin [43[ . Taking into account 
the non-flow contribution, the numerator of Eq. (|5J| can 
be written as follows [3(1 EH : 



(Q 2 ^>t)> = (^cos[2(^ T -4> k )\ 



(10) 



Mv 2 (pt)v2 



nonflow 



where 4> PT is the azimuthal angle of particles from a given 
Pt bin (v.2 i in Eq. ((5J) and the sum goes over all tracks 
k in an event used to determine the flow vector (Q2 in 
Eq. @). The angled brackets denote averaging over the 
events. The first term in the right-hand side of Eq. (fit))) 
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FIG. 4: (Color online) Charged hadron U2{FTPC} (closed circles) and v 2 {AA — pp, FTPC} (open circles) as function of pr in 
^jv" = 200 GeV Cu+Cu collisions for centrality bins: (a) 50 - 60%, (b) 40 - 50%, (c) 30 - 40%, (d) 20 - 30%, (e) 10 - 20% 
and (f) — 10%. The percentages refer to fraction of most central events. The error bars are shown only for the statistical 
uncertainties. 
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FIG. 5: Ratios of v 2 {AA - pp, FTPC}/v 2 {FTPC} for charged hadron as function of p T in ^/s^F = 200 GeV Cu+Cu collisions 
for centrality bins: (a) 50 - 60%, (b) 40 - 50%, (c) 30 - 40%, (d) 20 - 30%, (e) 10 - 20% and (f) - 10%. The percentages 
refer to fraction of most central events. 



represents the contribution from elliptic flow. v 2 (pt) is 
the value of elliptic flow at a given px- ¥2 is the ellip- 
tic flow on average for all particles used in the sum of 
Eq. (|10[) . The multiplicity of particles contributing to 
the sum is denoted by M. All other correlations subject 
to non-flow go to the second term in the right-hand side 
of Eq. (fit))) . It is assumed that the quantity {Q 2 u 2 Apt)} 
in p+p collisions can be used to estimate the non-flow in 
AA collisions MS- 



Mv 2 {p T )v 2 =(Q2U2i{pT)} 



AA 



(Q2U%i(pT))pp (11) 



Dividing both sides by 2y/ {Q2Q2*) aa as in Eq. flS) gives 



v 2 {AA - pp}(p T ) = 



(Q2ll2,i(PT))AA-{Q2U2 :i {pT)) r 



(12) 

because 2^J(Q a Q%*) A a = 2^(M/2)^(M/2)W = Mi*. 

Comparing p+p and AA collisions, one might expect 
some changes in particle correlations: there could be an 
increase in correlations due to a possible increase of jet 
multiplicities in AA collisions or, conversely, some de- 
crease due to the suppression of high px back-to-back cor- 
relations [44]] . On the other hand, AA collisions exhibit 
long 77 range correlations (the "ridge" ) [M, S] , which are 
not seen in p+p collisions and the origin of which is under 
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FIG. 6: (Color online) Charged hadron V2 integrated over pr and rj vs. centrality for the various methods described in the 
text in ^/snn = 200 GeV and 62.4 GeV Cu+Cu collisions. The error bars are shown only for the statistical uncertainties. 



TABLE I: Summary of systematic errors of v% due to the 
reconstruction procedure of strange hadrons in Cu+Cu col- 
lisions at t/snn = 200 GeV. 

K° s A 
Centrality Background Cut criteria Background Cut criteria 
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investigation [47[. Thus it is difficult to make an accu- 
rate estimate of non-flow contributions. The fact that at 
high pt (pt > 5 GeV/c) the p+p results are very close 
to central Au+Au [40, |43| suggests that the uncertain- 
ties are relatively small. In the following we estimate the 
systematic uncertainties arising from non-flow contribu- 
tions. We use v 2 {AA-pp, TPC} and v 2 {AA-pp, FTPC} 
to denote v 2 {AA - pp} calculated with TPC and FTPC 
flow vectors, respectively. 

III. SYSTEMATIC UNCERTAINTIES 

Non-flow is one of the largest uncertainties in elliptic 
flow measurements. As we mentioned in Section TlI C[ this 
effect can be investigated by comparing the azimuthal 
correlations measured in Cu+Cu collisions to those in 
p+p collisions. The event average of the sum of the cor- 
relations is given by Eq. (|TU1) . 

Figure [5] shows the azimuthal correlation, Eq. (ITTJl) . 
as a function of pt for the — 60% centrality range in 



Cu+Cu collisions at ^snn = 200 GeV, compared to p+p 
collisions. As we can see, the azimuthal correlations in 
Cu+Cu collisions, shown as solid squares, increase with 
Pt and then saturate above 2 GeV/c while those in p+p 
collisions, shown as open squares, monotonically increase 
with pt in the case of the TPC flow vector. With the flow 
vector determined from FTPC tracks the azimuthal cor- 
relations around midrapidity in p+p collisions are small 
when pt is less than 4 GeV/c. It means that one strongly 
reduces the non-flow effects with the FTPC flow vector 
relative to the one seen with the TPC flow vector. 

In order to illustrate the sensitivity to non-flow for the 
various flow analysis methods, we first analyzed h^ 1 el- 
liptic flow in the 60% most central Cu+Cu collisions 
at y/sNN = 200 GeV. As shown in Fig. [3] (a), the fact 
that t; 2 {TPC} is significantly larger than i> 2 {FTPC} in- 
dicates a larger non-flow effect in t> 2 {TPC}. With the 
large r\ gap between West and East FTPCs, non-flow ef- 
fects due to the short-range correlations arc reduced in 
u 2 {FTPC}. u 2 {FTPC} saturates sXp T ~ 2.5 GcV/c and 
then falls off slightly up to pt ~ 4 GcV/c. In order to 
estimate the remaining non-flow effects in v 2 {FTPC}, 
we subtract the azimuthal correlations of p+p collisions 
from those in Cu+Cu collisions according to Eq. (fl"2"|) . In 
Fig. [3] (a), v 2 {AA - pp,FTPC} is close to v 2 {FTPC} 
in the region pt < 4 GcV/c. To quantitatively il- 
lustrate non-flow systematic uncertainties, Fig. [3] (b) 
shows the ratios of v 2 {AA - pp,FTPC} to v 2 {FTPC}, 
v 2 {AA-pp, TPC} to v 2 {AA-pp, FTPC} and u 2 {FTPC} 
to w 2 {TPC} as a function of p T . w 2 {FTPC}/^ 2 {TPC} 
shows that non-flow in u 2 {TPC} increases from 20% at 
Pt ~ 0.8 GeV/c to 40% at p T ~ 3.5 GeV/c. Based on the 
comparison between v 2 {AA-pp, FTPC} and w 2 {FTPC}, 
the residual non-flow in u 2 {FTPC} is less than 10% below 
Pt ~ 4 GcV/c. We also checked the v 2 {AA — pp} calcu- 
lated with the TPC flow vector. Beyond pt ~ 3 GeV/c, 
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FIG. 7: (Color online) Charged hadron V2 as function of pt for 50 — 60% (solid circles), 40 — 50% (solid squares), 30 — 40% 
(solid triangles), 20 — 30% (open circles), 10 — 20% (open squares) and — 10% (open triangles) in ^/sjvjv = 200 GeV and 62.4 
GeV Cu+Cu collisions. The error bars are shown only for the statistical uncertainties. 
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FIG. 8: (Color online) (a) Charged hadron V2 as a function 
of pt in Cu+Cu collisions. The results from ^/sjvjv = 200 
GeV and 62.4 GeV are presented by open symbols and closed 
symbols, respectively, (b) Ratios of the V2(pr) from ^/snn = 
62.4 GeV to 200 GeV. The error bars are shown only for the 
statistical uncertainties. 



v 2 {AA—pp, TPC} seems systematically lower, but within 
errors it is similar to V2{AA — pp, FTPC}. This shows 
that most of the non-flow is eliminated by subtracting 
the azimuthal correlation in p+p collisions, validating our 



earlier assumption. 

To illustrate the ccntrality dependence of the system- 
atic uncertainties, Fig. @]shows V2{FTPC} and v 2 {AA — 
pp, FTPC} as a function of pr for six centrality bins. Ra- 
tios of v 2 {AA-pp, FTPC} to w 2 {FTPC} for each central- 
ity bin are shown in Fig. [5] from (a) the most peripheral 
bin 50-60% to (f) the most central bin 0-10%. For each 
centrality bin, the ratio falls off slightly as pr increases. 
For the two peripheral bins 50—60% and 40—50%, the ra- 
tios drop faster than in the other bins, indicating larger 
non-flow contributions in t>2 {FTPC} (pt) m peripheral 
Cu+Cu collisions. Figure [5] shows charged hadron v 2 in- 
tegrated over p T (0.15 < p T < 4 GeV/c) and 77 (\r)\ < 1.0) 
vs. centrality for the various methods. It is clear that 
V2{TPC} is much higher than for the other methods, es- 
pecially for the peripheral collisions. 

To summarize the non-flow systematics we employed 
the Scalar Product method with TPC and FTPC flow 
vectors for /i ± in Cu+Cu collisions at ^snn = 200 GeV. 
The results for the 60% most central events are shown in 
Fig. [3] «2{TPC} has large non-flow contributions while 
«2{FTPC} eliminates most of the non-flow. In what fol- 
lows, we will report our results in term of ^{FTPC}. 
For simplicity v 2 denotes u 2 {FTPC} except when the 
flow method is explicitly specified. With the assump- 
tion of pure non-flow effects in p+p collisions, we use 
v 2 {AA — pp, FTPC} to estimate non-flow systematic er- 
rors in v 2 {FTPC}. Ratios of v 2 {AA - pp,FTPC} to 
V2{FTPC} are shown for the 60% most central events 
in Fig. [3] (b) and six centrality bins in Fig. [S] The ra- 
tios show that non-flow effects increase with pr for all 
centrality bins and non-flow effects are larger in more 
peripheral bins. In order to estimate the non-flow sys- 
tematic error in V2{FTPC}, we fitted a constant to the 
ratio v 2 {AA~ pp,FTPC} /v 2 {FTPC} in the p T range (0, 
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FIG. 9: (Color online) V2 of K s (open circles), A ( solid squares), H (solid triangles) and (open stars) as a function of pr 
for (al) - 60%, (a2) - 20%, (a3) 20 - 60% and as a function ofm T -m for (bl) - 60%, (b2) - 20%, (b3) 20 - 60%. 
For comparisons, the results from ideal hydrodynamic calculations [49L l5(i| are also shown. At a given pr, from top to bottom, 
the curves represent n, K, p, cj>, A, 5 and Q. When pr is converted to hit — m, this mass hierarchy is reversed in the model 
results. All data are from tJsnn = 200 GeV Cu+Cu collisions. The error bars are shown only for the statistical uncertainties. 



4 GcV) . We take the numerical value of this constant as 
the estimate of the systematic uncertainty. The result- 
ing non-flow systematic error is minus 5% for — 10%, 
10 - 20%, 20 - 30% and 30 - 40% collisions and minus 
10% for 40 - 50% and 50 - 60% collisions. Although for 
Kg, 4>, A and S v 2 the non-flow effects may be different, 
as we don't have enough statistics to repeat the analysis, 
we simply assume a similar magnitude for the non-flow 
systematic error. 

The other systematic uncertainties in the v 2 analysis 
procedure are studied as follows. We estimate the sys- 
tematic errors from the shifting method for the FTPC 
event plane by comparing v 2 using different maximum 
harmonics in Eq. ([TJ and find the systematic errors are 
less than 1%. The systematic errors in K% and A v 2 re- 
sulting from the background uncertainty and topological 
cut criteria are estimated using the Event Plane method. 
The uncertainty due to the background subtraction is es- 
timated as the relative differences in v 2 from fitting the 
background using second and fourth order polynomials. 
The systematic uncertainty is also estimated by varying 
the cut parameters. The systematic errors for K s and 
A from the background uncertainty and the cut criteria 
are summarized in Table HI From Ref. [3l[, the estimated 
systematic uncertainty of A from feed-down is less than 



2% 



IV. RESULTS AND DISCUSSION 
A. Charged hadrons 

Flow results for charged hadrons were determined us- 
ing the Scalar Product method Eq.(8) with the flow vec- 
tor derived from the FTPC tracks. A comparison to 
v 2 {AA - pp, FTPC} Eq. ([12]) was used to estimate the 
systematic error. Figure [7] shows v 2 {pt) of h for six 
centrality bins from Cu+Cu collisions at ^/snn = 200 
and 62.4 GeV. For a given centrality bin, v 2 (j>t) initially 
increases with px- At higher px {pr > 2 GeV/c), v 2 
appears to saturate or decrease. v 2 {pt) m more periph- 
eral collisions increases faster and reaches higher values 
as expected for the larger eccentricity. 

At low pt, the increase of v 2 (pt) with px is consistent 
with predictions from ideal hydrodynamic calculations, 
which will be shown in Fig.[5]for identified particles. The 
model predicts that v 2 continues increasing beyond pr ~ 
2 GeV/c. The observed saturation or decrease of v 2 (px) 
indicates that the model is not valid in this region. One 
expects that the model should be valid up to higher px in 
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FIG. 10: (Color online) Same as Fig. [9] (al) and (bl), but expanded for the low px and tut — m regions. The data points with 
large errors have not been plotted. At a given pr, from top to bottom, the curves represent the ideal hydrodynamic calculations 
for 7r, K, p, <f>, A, a and f2 [49l. [50l] . When pt is converted to uit — m, this mass hierarchy is reversed in the model results. All 
data are from — 60% Cu+Cu collisions at ^snn — 200 GeV. The error bars are shown only for the statistical uncertainties. 



a system with larger densities and larger volumes. This 
was observed in 200 GeV Au+Au collisions [Hf where 
V2(pt) of saturated at higher px in more central col- 
lisions. However, we do not observe the strong centrality 
dependence of saturation pr for 200 GeV Cu+Cu colli- 
sions. 

Figure [8] shows the comparison of V2 for /i ± from ^/s^n 
= 62.4 and 200 GeV Cu+Cu collisions. The p T depen- 
dence of V2 at the two energies is similar. 



B. Identified hadrons 

The Event Plane method Eqs. ©-(O with event plane 
determined from the FTPC tracks was applied to Kg, A, 
S and (f>. The results are shown in Fig.[9]for the 60% most 
central events and also for the — 20% and 20 — 60% 
centrality bins at midrapidity \y\ < 1. Due to limited 
statistics, 3 and (j) are 0Iu y shown for the — 60% bin. 
The results from ideal hydrodynamic calculations (4j| [Ho| 
for each centrality bin are shown for 7r, K , p, 0, A, 3 and 
Q, which are displayed by the curves from top (bottom) 
to bottom (top) for the px (mx — m) dependence. 

The gross features of px dependence and hadron-type 
dependence are similar to those observed in 200 GeV 
Au+Au collisions [3l|. At low px, the hadron mass hier- 
archy (at a given px, the heavier the hadron, the smaller 
v 2(pt)) is reproduced by ideal hydrodynamic calcula- 
tions. (Sec Fig. [TU] (a) for clarity.) Multi strange-quark 
hadrons 5 and 0, which participate less in later hadronic 



interactions than do single strange-quark hadrons Kg 
and A, have sizable v-i- In particular, S is consistent with 
the mass ordering shown in Fig. |H](al). While the model 
can roughly reproduce the magnitude of the data for the 
60% most central events sample, there is an obvious dis- 
agreement in centrality selected data. The model under- 
predicts V2{px) in the — 20% bin while it over-predicts 
the data in the 20 — 60% bin. Effects not included in 
the model which may be relevant are geometrical fluctu- 
ations in the initial conditions and finite viscosity effects. 
It is unclear whether these effects can account for the 
difference between the model and data. 

At higher px, the hydrodynamic model mass order- 
ing breaks. V2{px) appears to depend on hadron type: 
V2(pt) is grouped by mesons and baryons, with mag- 
nitude depending on the number of quarks within the 
mesons or baryons. Over the entire px region, both the 
data and the model exhibit the same qualitative central- 
ity dependence as observed for 200 GeV Au+Au colli- 
sions [31j : the more peripheral the collision, the larger the 
V2 values. Compared to the results for 200 GeV Au+Au 
collisions 3l|, the splitting of Kg and A V2(px) is smaller 
in both the mass ordering region and the hadron-type de- 
pendence region. This indicates smaller collective flow in 
Cu+Cu than Au+Au collisions, which will be seen more 
clearly in Section llV Dl 

The transverse kinetic energy scaling first observed in 
Au+Au collisions is also tested in Fig. |H1 The results in 
Fig. [5] (al)-(a3) are re-plotted as a function of the trans- 
verse kinetic energy mx — m in Fig. |H] (bl)-(b3). The 



13 



0.08 
0.06 
0.04 
0.02 




t 1 1 1 1 r 



(a1) 0-60% 



■ n q scaling 



oK" s ■ A AS 



(b1) 0-60% 



0.5 1 1.5 



2.50 



0.5 



1 



P T /n q (GeV/c) 



1.5 2 
(m T - m)/n q (GeV/c 2 ) 



0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




I 1 1 r 

(a2) 0-20% 




t 1 1 r 

(a3) 20-60% 



/ 



t 1 r 

(b2) 0-20% 



Mil" 



J I I I UU I I I U I I L 



~i 1 1 r 

(b3) 20-60% 



m 

/ 



0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 0.5 1 1.5 

p^nq (GeV/c) (m T - m)/ n q (GeV/c 2 ) 

FIG. 11: (Color online) V2 /n q versus pr/riq (panels (al)-(a3)) and (my — m)/n q (panels (bl)-(b3)), where n q is the number 
of constituent quarks in the hadron. The parametrization Eq. ()13|) fitted to the data is shown as the dashed curves. All data 
are from ^/smn = 200 GeV Cu+Cu collisions. The error bars are shown only for the statistical uncertainties. 



quantity m denotes the rest mass of a given hadron. In 
the low rax — m region, v-i (jtit — m) is a linearly increas- 
ing function and independent of hadron mass. Transverse 
kinetic energy scaling holds in the region mj- — m < 0.8 
GcV/c 2 , as observed in Au+Au collisions jU [5l[. Cal- 
culations using ideal hydrodynamics arc shown in each 
panel as a function of rrvr — m. Contrary to the mass or- 
dering as a function of pr, the model shows the reversed 
mass ordering as a function of tjit — m: the heavier the 
hadron, the larger the vi{vtit — m) value. The results of 
Kg, A and !B exhibit tht — m scaling in each centrality 
bin, while the model does not show any scaling. Since no 
pion results are available, the scaling test of the data is 
not conclusive. All these effects can be seen more clearly 
in Fig. [101 

C. Quark- number scaling 

In Au+Au collisions in the intermediate px region, 2 
< pt < 4 GeV/c, the baryon-meson grouping of V2{pt) 
follows the number-of-quark scaling: the v-i of all hadrons 
fall onto a universal curve once Vi and pt are divide d by 
the number of quarks, n q , in a given hadron (3ll . |52| . 
The observed scaling can be explained by the coales- 
cence or recombination models [22l . l53l . [54|, indicating 



the constituent quark degree of freedom has been man- 
ifested before hadronization takes place. n q scaling is 
tested for various centrality bins in 200 GeV Cu+Cu col- 
lisions: V2(pt) and v-2.(jtlt — fn) scaled by n q are shown 
in Fig. QT] (al)-(a3) and (bl)-(b3), respectively. The rig- 
scaling formula of Ref. [5 51 ] , which can be written as 

^ iM = r°/ m / ~d (13) 

n q 1 + e -(Pr/«,-6)/c ^ I 

has been fitted to the data both in p? and tot — rn for 
each centrality bin. In this formula, a, b, c and d are fit 
parameters, n q is the number of quarks. The n q scaling 
is observed for pr/riq > 0.8 GeV/c, whereas it is seen for 
the entire (mT~m)/n q region. Below (mT~m)/n q ~ 0.4 
GeV/c 2 , the — m scaling which was established at 
low pt, now scaled by n q leads to the combined (mj — 
m) / n q scaling. The universal n q scaling of vi suggests the 
manifestation of early partonic dynamics in both Au+Au 
and Cu+Cu collisions. 



D. Centrality and system-size dependence 

The centrality and system-size dependence of vi is re- 
lated to the physics of the system created in high energy 
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TABLE II: Participant eccentricity e par t{2} and number of participants Apart from the Monte Carlo Glauber model |23, H3] 
and Color Glass Condensate (CGC) model [58M6ll ] calculations in Au+Au and Cu+Cu collisions at y / sjvw = 200 GeV. The 
quoted errors are total statistical and systematic uncertainties added in quadrature. 



Au+Au 



Cu+Cu 



Centrality e part {2}(CGC) 


Epart {2} (Glauber) 


Apart 





- 80% 


0.338 ± 0.002 


0.302 ± 0.004 


126 ± 8 





- 10% 


0.148 ± 0.001 


0.123 + 0.003 


326 + 6 


10 


— 40% 


0.353 ± 0.001 


0.296 ± 0.009 


173 ± 10 


40 


— 80% 


0.554 ± 0.002 


0.533 ± 0.018 


42 + 7 





- 60% 


0.336 ± 0.009 


0.350 ± 0.008 


51 + 2 





-20% 


0.230 + 0.010 


0.235 ± 0.008 


87 + 2 


20 


-60% 


0.434 ± 0.003 


0.468 + 0.016 


34+1 





-10% 


0.187 + 0.002 


0.197 + 0.002 


99 + 2 


10 


-20% 


0.281 + 0.002 


0.279 ± 0.008 


75 + 2 


20 


-30% 


0.360 ± 0.003 


0.369 ± 0.009 


54+1 


30 


-40% 


0.428 ± 0.002 


0.458 + 0.017 


38+1 


40 


-50% 


0.490 ± 0.002 


0.550 + 0.021 


26+1 


50 


-60% 


0.555 + 0.004 


0.643 ± 0.031 
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FIG. 12: (Color online) U2 scaled by participant eccentricity as a function of pt in y^ijvjv = 200 and 62.4 GeV Cu+Cu collisions. 
The error bars are shown only for the statistical uncertainties. 



nuclear collisions. In the ideal hydrodynamic limit the 
centrality dependence of elliptic flow is mostly defined 
by the elliptic anisotropy of the overlapping region of the 
colliding nuclei, and in the low-density limit by the prod- 
uct of the elliptic anisotropy and the multiplicity. Thus, 
the centrality and system-size dependence of elliptic flow 



should be a good indicator of the degree of equilibration 
reached in the reaction 28]. 

For a study of the centrality dependence of V2{jpt) in 
Cu+Cu collisions together with Au+Au collisions, we di- 
vide V2(j>t) by the initial spatial anisotropy, eccentricity, 
to remove this geometric effect. The participant eccen- 
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FIG. 13: (Color online) Centrality dependence of V2 scaled by number of quarks and participant eccentricity (t>2/(n g e par t{2})) 
for K% (left) and A (right) as a function of (mr — m) /n q in 0—10%, 10 — 40% and 40 — 80% Au+Au collisions (open symbols) [3l| 
and — 20% and 20 — 60% Cu+Cu collisions (solid symbols) at ^/sjvjv = 200 GeV. Curves are the results of n^-scaling fits from 
Eq. (| 13p normalized by £ par t{2} to combined Kg and A for five centrality bins. At a given pr, from top to bottom, the curves 
show a decreasing trend as A part decreases. The error bars are shown only for the statistical uncertainties. 



tricity is the initial configuration space eccentricity of the 
participants which is defined by |56f 



-part 



In this formula, al 



4«) 



^7, 



(14) 



(x)\ a\ 



(y 2 ) - (v) 2 and 

&xy = (xy) — (x)(y), with x, y being the position of the 
participating nucleons in the transverse plane. The root 
mean square of the participant eccentricity 



-part 



{2} 



-part/ 



(15) 



is calculated from the Monte Carlo Glauber model I29l . 
\El\ and Color Glass Condensate (CGC) model [58l - l6l| . 
(See TableUfor e pa rt{2}.) Since the FTPC event plane 
is constructed from the hadrons which have their origin 
in participant nucleons and the FTPC event plane resolu- 
tion is less than 0.2, what we actually measure is the root 
mean square of v 2 with respect to the participant plane 
[62^ . In this case, e pa rt{2} is the appropriate measure 



of the initial geometric anisotropy taking the event-by- 
event fluctuations into account j62H64| . Figure [T2l shows 
the centrality dependence of ^2(PT)/^part{2} for h in 
200 and 62.4 GeV Cu+Cu collisions. For a given central- 
ity bin, f2(pr)/£part{2} initially increases with px and 
then flattens or falls off at higher pr- After the geomet- 
ric effect is removed, the ordering of the distributions as a 
function of centrality, observed in Fig. [71 is reversed: the 
more central the collision, the higher the W2(PT)/e P art{2}. 
This suggests that the strength of collective motion is 
larger in more central collisions. 

To further study the centrality dependence of strange 
hadron v 2 , we normalized the n g -scaled values by 
£part{2} and plotted them as a function of (mj — m)/n q . 
The centrality dependence of Kg and A results are shown 
in Fig. [T31 The full symbols show from top to bot- 
tom the results from — 20% and 20 — 60% centrality 
Cu+Cu collisions. For comparison, the results from 200 
GeV Au+Au collisions HH are shown by open symbols 
in Fig. [T31 The results in Au+Au collisions are slightly 
different (~ 10% larger) from the previous published re- 
sults [3l|, which were calculated directly from the wide 
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FIG. 14: (Color online) Number of quarks and participant eccentricity scaled i>2 (V2 / (n q e pa , T t{2})) of identified particles as a 
function of (mr — m)/n q in — 80% Au+Au collisions (open symbols) [3l|] and — 60% Cu+Cu collisions (closed symbols) at 
yfsNN = 200 GeV. Circles, squares and triangles represent the data for Kg, A and S, respectively. The error bars are shown 
only for the statistical uncertainties. 



ccntrality bins. From top to bottom, the results are from 
- 10%, 10 - 40% and 40 - 80% centrality bins. For 
clarity, Kg and A results are shown in different panels. 
Curves represent n g -scaling fits from Eq. (|13[) normalized 
by £ P art{2} to the combined data of Kg and A for five 
centrality bins. For a given centrality, Kg and A results 
follow a universal curve, which means partonic collective 
flow is explicitly seen in the measured scaling with n q 
and £ P art{2}. For a given collision system, the stronger 
partonic collective flow is apparent as higher scaled Vi 
value in more central collisions. To study the system- 
size dependence of the scaling properties, the results from 
- 60% centrality Cu+Cu and - 80% Au+Au collisions 
are shown in Fig. 1141 The stronger collective motion in 
Au+Au compared to Cu+Cu collisions becomes obvious 
although the constituent quark degrees of freedom have 
been taken into account in both systems. 

In the ideal hydrodynamic limit where dynamic ther- 
malization is reached, the mean free path is much less 
than the geometric size of the system. The geomet- 
ric size of the system and the centrality dependence of 
flow is totally governed by the initial geometry (eccen- 
tricity) [28| . As there is no universal scaling with the 
eccentricity among either different collision centralities 
or different collision system sizes, this indicates that the 
ideal hydrodynamic limit is not reached in Cu+Cu colli- 
sions, presumably because the assumption of thermaliza- 
tion is not attained. In addition, W2/(^ij£part{2}) shows 
an increasing trend as a function of N part (See Fig. IT5|) . 
Table [TT| lists the values of eccentricity and N palt for the 



used centrality bins in Au+Au and Cu+Cu collisions. 
This suggests that the measured V2 is not only depen- 
dent on the initial geometry, but also on A^ part . 

Theoretical analyses found that the ccntrality and 
system-size dependence of vi can be described by a sim- 
ple model based on eccentricity scaling and incomplete 
thermalization. Within these models the lack of perfect 
equilibration allows for estimates of the effective parton 
cross section in the quark-gluon plasma and of the shear 
viscosity to entropy density ratio (t]/s) [33l. |65|). Thus, 
the V2 results from Cu+Cu collisions reported in this pa- 
per should allow extraction of 77/s and extrapolation to 
the ideal hydrodynamic limit. 



V. SUMMARY 

We present STAR results on midrapidity elliptic flow 
vi for charged hadrons /i ± and strangeness containing 
hadrons Kg, A, S and (f> from Cu+Cu collisions at ^/snn 
= 62 A and 200 GcV at RHIC. The centrality dependence 
of V2 for different system-sizes as a function of the trans- 
verse momentum pt is presented. To estimate the sys- 
tematic uncertainties, we studied various measurement 
methods. Below px ~ 4 GeV/c, non-flow correlations are 
reduced with the event plane constructed from hadrons 
produced in the forward regions (2.5 < \r]\ < 4.0). We 
obtained an estimate of the systematic uncertainties due 
to remaining non-flow contributions based on correlations 
measured in p+p collisions. 
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For a given centrality bin, px and hadron-type de- 
pendences of strange hadron i> 2 are similar to those 
found in Au+Au collisions [3l|: (i) in the low px re- 
gion, pt < 2 GeV/c, the hadron mass hierarchy is ob- 
served as expected in ideal hydrodynamic calculations: 
at fixed pr, the larger the hadron mass, the smaller the 
V2- (h) in the intermediate px region, 2 < pt < 4 GeV/c, 
i)2 as a function of either px or nix — m follows a scal- 
ing with the number of constituent quarks n q . Larger 
^2/(^g£part{2}) values are seen in more central collisions, 
indicating stronger collective flow developed in more cen- 
tral collisions. The comparison with Au+Au collisions 
which go further in density shows eccentricity scaled V2 
values depend on the system size (iV P axt)- This sug- 
gests that the ideal hydrodynamic limit is not reached in 
Cu+Cu collisions, presumably because the assumption of 
thcrmalization is not attained. 
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